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Abstract 

We develop a framework and unified, compact notation to include neutrino 
direct interaction effects with neutrino oscillations for a wide class of low- 
energy, effective four- Fermi interactions. Modest flavor- violating interaction 
strengths can make significant changes in the boundaries in AM^ sin^26 
and in tan^Ois ^ tan^d23 parameter space. We draw examples from the L — R 
symmetric and MSSM + ^ models and find the recently reported decay-at-rest 
jjL e transition probability can be described without oscillations m. L — R 
symmetric models. 

PACS: 13.15.+g, 14.60.Pq, 12.60.-i, 12.60.Jv. Keywords: neutrino, oscilla- 
tion; neutrino, flavor; muon-|-, leptonic decay; pi-|-, leptonic decay; interaction, 
neutrino p; numerical calculations, interpretation of experiments. 

1 Introduction 

Reconciling the solar atmospheric Q, and accelerator ^ || neutrino data 
is difficult within a purely 3-family oscillation description, which provides for 
two mass-squared differences, 3 mixing angles and a CP violating phase @, 0, §|. 
With this in mind we develop a unified framework to combine the oscillation new 
physics with, presumably related, direct-interaction new physics. We apply the 
general formalism to a number of cases with an eye to accommodating all of the 
positive signals for new phenomena within a three-family framework. We focus 
mainly on the picture where the favored |^ Ve ^ MSW||l^ oscillation is the 
solution to the solar neutrino observations, a <-> z^,- oscillation is the solution 
to the atmospheric neutrino anomaly and a direct interaction of the type 
Vp, + n ^ +p at the detector and/or 7r+ fi^ + Ve or /x+ e+ + i?e + z/^ at the 
source contribute substantially to the LSND signal[|[ Q for the 'Ve appearance" 
at the few tenths of a percent level. The direct interaction effects are flight- 
path independent and, moreover, they are clearly too small to explain the large 
solar and atmospheric anomalies that have been reported. At the same time, 
the strictly two-flavor t'/x oscillations with currently preferred[9] mixing 

and mass values, are negligible in the short-baseline experiments of the LSND 

^Email address: johnson@kuphsx.phsx.ukans.odu 
^Email address: mckay@kuphsx.phsx.ukans.edu 



1 



type. To complement the "two two-flavor", Ue ^ Vfi plus <-> Vr, picture 
just mentioned, we include a brief discussion of a simple "one-scale" version of 
three-flavor mixing within our framework |Q. We show the addition of modest 
direct interaction effects extends the allowed parameter space in an interesting 
way. 

A clean and convenient parameterization of the situation can be achieved 
by separating the problem into the production (source), propagation (with or 
without matter enhancement effects) and observation (detector) components 

First we introduce a working phenomenological model which includes 
lepton- flavor violating, four- Fermi interactions of purely leptonic and of semilep- 
tonic types.0 Next we present the setup to describe the production, propagation 
and detection and define the "appearance probability" for the situation where 
the neutrinos are relativistic, the source is a muon or charged pion, the prop- 
agation is in vacuum and the detector is a neutron or proton (or nucleus) tar- 
get. We then survey serveral patterns of dominant interactions suggested by the 
data. We compare these patterns with those presented by the minimal left-right 



symmetric model f5, ^ and the minimal supersymmetric model augmented 
by a minimal set of discrete R-parity violating terms [T^, p^ . Considering the 
LSND|, I result by itself, we make a study of the shift in allowed mass and 
mixing angle parameter space that is afforded by including direct interactions. 
In our concluding discussion, we point out several open possibilities. 



2 A Working Model 

To treat the bulk of the current accelerator data and to illustrate the general 
ideas, it suffices to take a model with two lepton families and one quark family 
with charged current, effective four-Fermi interactions. The model semileptonic 
(S) and leptonic (L) effective Lagrangian takes the form ig// = + L^, where^ 

2''/^GFOA{K'Xj,rAPhUja'ya) + h.c. (1) 

= 2^/^GFF^^^^J,rAPhUj,MikrAPhamb'yb)\ (2) 

where the coefficients K and F give the amount of admixture of /i — L, R 
projection, the type of Dirac matrix structure and the lepton flavors contained 
in Leff. Possible crossed terms between L and R give contributions to the 
unpolarized cross sections that are proportional to m^/'s and we ignore these 
effects here. The hadronic operators are of form Oa = dTAiaPL + bPFi)u. All 
repeated indices are summed: A — S,V,T (matching Dirac gamma matrices 

^The questions of the detailed form and origin of the neutrino mass matrix are not ad- 
dressed here. 

*The factor of ^ in the definitions Pj^ = ^ (1 — 75) and Pji = ^ (1 -I- 75) leads to the factor 

2^/2 ^]-ign combined with the conventional factor in the four-Fermi interactions. Neutral 

V 2 

current SM terms are also contained the term by Fierz symmetry of {V — A) X (V — A). 
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Table 1: Some relevant expressions for K and F coefficients in the SM, L — R 
and MSSM + ^ models. The L-RanA^ entries should be divided by Gi^23/2. 











^yii22 


Fv211m 


SM 


1 








1 





L-R 











/l2l/li2 


/ill /l2m. 


n 




'^21fc-^llfe 


/ ^ 

'^311-^231 


-|Al23|^ 





^ M?, 

k d 




<^3 


R 



1, 7 and cr), h — L, R, lepton flavor indices i, j, fc, m and mass eigenstate indices 
a, 6. The C/ja are the unitary transformation matrices between the flavor states 
and the mass eigenstates. 

In a two-flavor mixing scheme where I'e and stand for the fields whose 
Pl projections couple with e and fi to the W-boson, we express i^e and in 
terms of the "mass eigenfields" vi and 1/2, and write z/g = cosO vi + sinO vi 
and fc'p — —sin9 vi + cosO 1^2 ■ Several examples of the coefficients K and 
F are tabulated for the SM and two familiar extensions in Table (|l]). In 
the ultrarelativistic limit, an approximate Fock space of flavor states can be 
constructed |l3|, while the > and |z^2 > states propagate with the simple 
plane wave factors exp i(x • p — Ei_2t), leading at order (Am^)/i? to the usual 
mass, energy and distance dependent argument of the oscillation probability. 
The task at hand is to fold this oscillation effect with the source and detector 
lepton-flavor violating effects to produce an "appearance probability" in terms 
of a convenient, uniform parameterization. 



3 Analysis of an Accelerator "Appearance" Pro- 
cess 

Consider 7r+ — > i/ + /x"*" decay to be the source of neutrinos, where the transition 
matrix element to a mass eigenstate \vb > reads M^ —< i^f,, ^+|L'^|7r+ >, with 
L'^ given in Eq. (|l]). Consider the reaction z/ + 7Vi(>l, Z- 1) ^ iV/(A, Z) + to 
be the detection mechanism, with matrix element M^ =< e^, Nf\L^\Ni, Vb >. 
We may now write the /i — > e transition probability in the ultrarelativistic 
neutrino limit as|p|, P^^e ^ |S6Mj^e"'-^''*M^fp, where the sum is over 
mass eigenstates and the equality holds to leading order in m^/E. Spelling out 
this probability in the notation of Eq. (|]) , we have the expression 

p^^e ~ {2y'GFf\j:t{K'x^,rK^,kif^rAPkU,t>^brx 

eVsPhUkbVbe-''^^^ < 0\Oa\tt+ >< Nf\OB\Ni > p. ^ > 

The sums over repeated indices A, B, h, j, k and a sum over spins are implicit. 
In the same fashion as above, the /i+ e'^ + v + 9 process can be considered to 
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be the source of neutrinos, with the effective Lagrangian replacing in the 
source matrix element. The "appearance" sought is then v + Ni{A^ Z + 1) ^ 
Nf{A, Z) + e+. The counterpart to Eq. (||) reads: 

{erAPhU„,,i^,re-'^^\erBPhUkd'^dr < Nf\OB\m > |2. (4) 

The expression in Eq. is appropriate to the LSND decay-in-flight (DIF) 

study and Eq. to their decay-at-rest (DAR) study with Ni = p and Nf = n. 
To make the notation more concrete, let us illustrate it with several examples. 

4 Purely Leptonic New Physics 

Our first example is applicable to the DAR situation, where fj,'^ decay provides 
the D^s at the detector. We assume only the lepton source contains new, lepton- 
number violating interactions, and that they are represented in the effective 
Lagrangian, Eq. (||), by {V — A) x {V A) terms. The lepton-flavor violating 
case of interest to us has i — 2, j = f, k = 1 and m = f , 2 or 3. Applying 
Eq. (^) to this situation, one reads off, for ^+ e^vv and i?e + ^{Z) 
e+ +N{Z~1), 

c d 

(5) 

With the new physics confined to the lepton-number violation in the decay 
as described, we have Ky-^^ = 1 = Py22ii ^'^^ ^v'2iim ^- The factor in braces 
is exactly the one appropriate to the unsupressed sequence /^^ e^v^Ve and 
i?e + N{Z) N{Z — f ) -I- e+. The P^j_,e probability can now be read off from 
Eq. (^), and for two flavor mixing we obtain 

Pfj^—te — E [^^^^^V'm + (1 ~ ta'n?'ip,ji)si'n?29sin'^x] 

m—e,fi,T . 

[ta,n ipi,sin4:9sin'^x ^ ) 

-\- 4:tamp^sin26sinipi,sinx{cosipeCosx — cos29 sirup esinx)], 

where 2x — {E2 — Ei)t and where Pv'2iim = tcLmpraR^^'^"^ has been defined so 
that a uniform "all angles" representation of the appearance probability could 
be achieved. The limit — *■ 0, keeping only an m = /i value and keeping 
leading powers of V'm << Ij leads to the formula of Ref. |0 for the case that 
the direct lepton-flavor violation is only in the /i-decay source interaction. 

In Fig. @) we show the P^^e = 0.003 ± 0.001 1|| contours in the sin^2e, 
AM^g parameter plane for V'e = Pv'2111 values of and 0.04, where — 'tpr = 
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Figure 1: The LSND la allowed region (P = 0.002-0.004 for DAR) for no direct 
interaction new physics (solid curves) and purely leptonic new physics strength 
tjj = 0.04 (dashed curves). Nonzero new physics phase makes the curves slightly 
lower for AM^ < 1.8eV^. Also included are limits from Bugey (dash-dotted 
curve) and E776 (dotted curve). 

Table 2: Relevant new physics couplings and upper bounds. Sources of the 
bounds are other than neutrino oscillation experiments. 
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ipm = is assumed.^ The minimal L — R model [|T5[ is relevant to the case 
under consideration, since it generates lepton-flavor violating, purely leptonic, 
four-Fermi, {V — A) x {V — A) interactions by exchange of an iso-triplet Higgs 
particle. There is no {V-A)x [v-A) fi+ e++De+'^x effect in the MS'S'M-h^ 
model, contrary to the assertion in Ref. [14]. In the notation of Eq. (|^) we 

have 2'^/^GFFy^j^^ = ^Mi'"" ^ tti^ identification of the relevant coefficients in 

the effective Lagrangian, Eq. (||). Yukawa coupling matrices are designated by 
hij = \{hij + hji) and Ma^^ is the mass of the singly-charged member of the 
Higgs triplet, following the notation of Ref. The ranges of allowed values 

of for the L - i? model are shown in Table (||). It is surprising that 

a considerable flexibility in choosing relevant L — R model parameters is still 
allowed by all experimental constraints. 

5 The Same New Physics at Source and Detec- 
tor 

A case that injects the same new semileptonic physics at source and detector 
in the DIE situation is exemplified by the pure {V — A) x {V — A) choices 
Ky22 — Kyii — 1' which follow from the standard model, and Ky^^ — Ky^2 = 

^ Since we are showing only t.re nHs and not fits, we adopt tiie Gaussian "toy model" used 
for illustrating oscillations in Ref. |20{ with bo and o-j, appropriate to LSND. 
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Figure 2: The LSND la allowed region (P = 0.0015 - 0.0037 for DIF) for no 
new direct interaction physics (solid curves) and for P — 0.0026 with = 0.01 
(dotted curve), ip = 0.02 (dot-dashed curve) and -0 = 0.03 (dashed curve) when 
there is the same new physics at the source and detector. 

tantpe^^'^ which defines the new physics.^ The Pe^fj. expression that results 
when two-flavor mixing is assumed reads 

Pfj,^e{DIF) = Atari^^ cos^x + sin^29sin^x — 2tan'4'sin29sin2ipsin2x (7) 

for DIF]^ On the other hand, 

P^^e{DAR) = tan'^tp + sin'^26sin'^x {1 - tan'^ip) . 

+ tan^sin29 [—sin2(psin2x + 2cos2(fcos26siii?x] 

for DAR, where only the detection is aff'ected by the new interaction. In writing 
Eq. (^, we have assumed that there is no lepton-flavor violating new physics 
interaction in /i-dccay. If sin^x « sin^ip^ then P^^e{DIF) ~ AP^^e{DAR), 
so a combination of and effects would have to be invoked to produce the 
P^^e{DIF) ~ P^,^e{DAR) ~ 0.003 resuh reported by the two LSND experi- 
ments. Fig. d) shows plots of the P^,^^{DIF) = 0.0026 ± 0.0011 boundaries 
in the AM^^ ^ sin^29 parameter plane for several choices of ip with (p = 0. 
For smaller ^ values, ip <, 0.02, the modifications to the ip = curve have the 
same general character as those of Fig. (P, moving the curves leftward and 
downward. A dramatic change occurs when ^ 0.03, however. In Eq. (^ the 
large term Atan^ipcos^x tends to "overshoot" the input value of Pfj,^e as tan^ip 
grows, so X needs to increase to compensate this growth. But the second term 
grows as x grows, so sin^26 must control this growth and therefore a maximum 
sin^26 value is established for larger tari^ip values. It is the factor four in the 
first term of Eq. (^, reflecting the conipoimd effects of source and detector, 
that produces the new effect at rather modest tan^ip values. 

To make the connection to the minimal ft SUSY model, we read off 

Ky2i — Kvi2 — Yl, ^^M^^^^'' ' 23/2G fro™ the semileptonic effective Lagrangian 

that follows from exchange of the c?^ squark, with k the family index, in the 
notation of Ref. |]l9j. In this model, the neutral current e ^ /i transition 
operator has the same coefficient, and the bound from /iTi eTi ||2l|, ^ is 
severe, as listed in Table (H), making this case academic in the minimal ft 
SUSY model. At the moment we cannot offer a well-motivated model with the 

^We use the same symbols ip and ifi to parameterize the different cases for economy of 
notation. The context makes clear that these are independent parameters for the different 
cases. 

^On this point, we completely disagree with Ref. |l4], which asserts that there is no effect 
when the physics at source and detector are the same. 
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structure Ky^i = ^vi2' but we present the above to illustrate the structure of 
the case where the source and detector new physics is the same. 

6 New Physics Only at the Detector 

Because the DAR and DIF experiments of the LSND collaboration have different 
backgrounds and different systematic errors, the reported agreement between 
the results of the two experiments suggests a common origin[||. We point out 
here that a lepton-number violating effective interaction that is only operative 
in the detection processes tends to produce this effect. 

Consider the effective Lagrangian case with ify22 = ^Lii = 1 the stan- 
dard model, Ky^2 7^ ^ tanipe'^^'^ and all other coefficients negligible for the 
new physics. We find that, in the two-neutrino mixing case, the same param- 
eterization, Eq. (||), applies to both the DAR and DIF cases even with the 
direct, lepton-number violating interactions effects combined with those of the 
oscillations . In principle the agreement between the DAR and DIF measured 
oscillation probabilities could be represented by tan'^il', a purely direct 
interaction effect with tan'^ijj ~ 0.003 

In surveying our example models, the minimal SUl{2) x SUjf {2) x U{1)b-l 
weak interaction model and the MSSM without i?-parity, one finds that the 
above parameter choice Ky-^^^^ and all other new physics K^-j = 0, does not 
occur. The instructive case 0, ^^521 — and with Ky^^ = by choice, 

does arise from the i?-parity violating interactions, however, and we summarize 
this case next. 

7 Different New Physics Lorentz Structures 

Up to this point, we have considered only cases where the new interaction, 
charged current effective four-Fermi Lagrangian has the same {V — A)y. {V ~ A) 
structure as the standard model. This ensured all matrix elements of interest 
had the same structure and the dynamics of the source and detector physics 
could be lumped into an overall factor that is normalized away in the definition of 
the lepton "flavor change" , or "oscillation" probability. If the Lorentz structure 
of the charged current x charged current interaction is not the same as that of 
the standard model, then the analysis of P^^e necessarily brings in details of 
the matrix elements at the source and/or detector. 

In the situation mentioned above where ^^512 ^ is the only relevant new 
physics, and the scalar interaction is of the form {S — P)iepton x {S + P)quark, the 
matrix elements of the new physics and those of the standard model operators 
are distinctly different. Dividing out a common factor from the source, the DAR 
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case looks schematically like, in the two-flavor mixing case, 



Pt.-.e ~ \nxPL'^<Nf\0^\N,>{~stn2esinxe") 

+ K^^^ePLU < Nf\Os\Ni> {\- 2cos^esinxe")\^ . ^ ' 

The new physics, lepton-number violating interaction effects depend both on the 
parameters (^^512, S and x) and on the comparative role of the matrix elements: 
e^fxPhV < Nf\0^\Ni > in the standard model term and cPli^ < Nf\Os\Ni > 
in the new interaction. With Os = uPud, as in the i?-parity violating model 
where K^^^ = X',\,X23i/{MIGf2V2),K^^, = and K{^,, = K^\^ « K^,„ 
we can evaluate the spin-averaged transition probability for ai^ + n— > e + p 
detector transition to illustrate the DAR application. Simplifying the /?-decay 
matrix element to pure y — ^ for illustration purposes, we find 



Pf,^e sin^ 2esin'^x + ji^liaP^T^ — " sin^2esin'^x), (10) 



where Mn is the mass of the nucleon, s — (P^ + P„)^ and t — (P^ — Pe)^- 
In writing Eq. (|l^) we have dropped the cross term, which is proportional 
to TOe, and we have not displayed a final state phase space integration. The 
matrix element effects are clearly crucial, as is shown by the suppression of the 
cross term and the nontrivial nature of the relative kinematical factor between 
the {V - A) X {V - A) standard model structure and the {S - P) x {S + P) 
new physics structure. For energies appropriate to LSND, the kinematical fac- 
tor ratio that multiplies the new physics term is small (of order |</m^| ~ 
(SOMeV/lGeV)'^), so it is unlikely that such terms can contribute significantly 
to a Pfi^e signal. For completeness, we mention that the MSSM with ft, where 
l'^3ii'^23i Gf2\/2 < 2 X 10^^|22 , requires a kinematical enhancement to 

^ '3 



bring it up to the order of 10~ relevant to the LSND reported result. The 
important lesson to stress, however, is that one must look at the detailed ma- 
trix element evaluations to draw conclusions in these cases where the Lorentz 
structure of the effective Lagrangian's new physics terms is different from that 
of the standard model. 



8 Three— Family Mixing Application 

Our illustrations so far have relied on two-family mixing parameterizations, 
which are appropriate to the M3 >> M2 >> Mi type of hierarchy of neutrino 
masses. The framework we present applies to any number of families, so we 
discuss a simple, "almost viable" version of three-family mixing as our last 
illustration |Q. The starting, simplifying assumption is that Mi ^ M2 << 
M3, which reduces the three-family mixing to a "one-scale" problem; namely 
AMfg ~ AMI3 = AM^ with AMf^ ^ 0. Returning to the DAR case with new 
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Figure 3: For three-flavor, one-scale mixing, the LSND Icr allowed region (DAR) 
for no new direct interaction physics (solid curves) and purely leptonic new 
physics strength — 0.02 (dotted curves) and ip = 0.04 (dashed curves). 

physics at the source and a restriction to m = 1 in the flavor sum to reduce 
clutter in the notation, we have 

Pe^ = I - 2iU;3Ui3sinxe-'^ + tan^e'^^'^{l - 2i\Uis.\^ sinxe-''')\^ , (H) 

where JV2iii = tanipe'^'"^. With the parameterization U23 ~ sin923Cos9i3 and 
U13 — sin9i3e~^^^^ [pO| , one can write out an expression similar to, but even 
more obscure than Eq. (H). The form shown in Eq. ( pi]) makes the point 
that the first term alone gives the result of Ref. |^ , and the new interactions 
produce a leading tan^ijj term plus other terms similar to those in Eq. (||). The 
flexibility afforded by the inclusion of the direct interaction is shown in Fig. 
(||), where the effect of increasing ip on the tan^9i3 vs. tan^923 plot is shown to 
be significant for i/j ~ 0.04. 

9 Discussion 

We have developed and applied a compact analysis that treats neutrino fiavor- 
changing phenomena - oscillations plus direct interactions - in a unified manner. 
We showed standard oscillation parameter plots can be shifted significantly for 
rather modest values of -0, an angle parameterizing direct interactions. We also 
showed the positive DAR result reported by LSND[^ can be accounted for by 
direct interaction effects in the L — R model. 

Open possibilities for further applications include the unified treatment of 
the neutrino masses and direct interactions within a model like the MSSM + ft 
model and the applications to general three-family and three-family plus sterile 
neutrino pictures. 

Whether direct lepton-fiavor changing effects turn out to be important or 
not, the framework offered here should provide an efficient means to survey the 
possibilities. 
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